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Ferromagnetic Polarons in Manganites 
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Using the Lanczos method in hnear chains we study the double exchange model in the low con- 
centration limit, including an antiferromagnetic super-exchange K. In the strong coupling limit we 
find that the ground state contains ferromagnetic polarons whose length is very sensitive to the 
value of K/t. We investigate the dispersion relation, the trapping by impurities, and the interaction 
between these polarons. As the overlap between polarons increases, by decreasing K/t, the effec- 
tive interaction between them changes from antiferromagnetic to ferromagnetic. The scaling to the 
thermodynamic limit suggests an attractive interaction in the strong coupling regime (J^ > t) and 
no binding in the weak limit (Jh — t). 



I. INTRODUCTION 

The discovery of "colossal" magnetoresistance (CMR) 
[01 together with its many unusual properties has re- 
ceived considerable attention lately. The materials that 
exhibit this phenomenon are under much experimen- 
tal investigation due to their technological applications. 
These perovskites are ferromagnetic oxides of the form 
Ai_j;B:rMn03 (where A=La,Pr,Nd; B=Sr,Ca,Ba,Pb) 
1^,^. Experiments have revealed very rich phase dia- 
grams interpreted in terms of ferromagnetic, antiferro- 
magnetic, canted, polaronic, and inhomogeneous phases. 
Charge ordered phases have also been found in these 
compounds. 

The phase diagram, as a function of concentration x, 
temperature, magnetic field, or magnitude of the relation 
between super-exchange and hopping interactions is not 
yet completely understood for the different compounds. 
A clear asymmetry is observed between hole (x < 0.5) 
and electron doped {x > 0.5) regimes. A metallic ferro- 
magnetic phase can be reached by hole doping of the par- 
ent compound LaMnOa, by substituting La for divalent 
alkalis, Pb, or by stoichiometry changes. On the other 
side, electron doped systems exhibit charge ordering or 
quasi-ordering and antiferromagnetism. 

In the hole doped perovskites, the strong correlation 
between itinerant carriers and localized spins gives rise 
to competing magnetic interactions. This competing in- 
teractions plus the effects of disorder and lattice degrees 
of freedom, could give rise to inhomogeneities of vary- 



ing scale: microscopic polarons, mesoscopic droplets or 
macroscopic phase separation. A number of experimen- 
tal data indicating the different scales in the hole doped 
regime can be found in Refs. |p^ , p| for polarons, |12| for 
droplets and Q for phase separation. 

Much less is known about the electron doped com- 
pounds. This lack of initial interest can be explained 
by the fact that manganites such as Bio.25Cao.75Mn03 
or Ao.25Cao.75Mn03 are characterized by charge order- 
ing phenomena and do not exhibit magnetoresistive ef- 
fects. However the recent discovery of new ferromag- 
netic electron doped manganites, Cai_j;Bj;Mn03 and 
Cai_j;Bi^Mn03 for x = 0.10 ]l^ , p^ , has increased the 
interest for those compounds since it opens the possibil- 
ity of magnetoresistance effects also in this concentration 
limit. Cao.9Bio.iMn03 and Cao.9Euo.iMn03 exhibit in- 
deed negative magnetoresistance but with a much smaller 
magnitude than hole doped manganites. Long range 
charge ordering has recently been observed with elec- 
tron diffraction in Lai^^^Ca^MnOa {x > 0.5) §-|,|l6|,|l3 . 
The magnetic susceptibility has a pronounced inflection 
at the charge ordering temperature, resembling that as- 
sociated with a conventional long range antiferromag- 
netic transition . Neutron scattering measurements 
in Bii_^Caa;Mn03 (0.74 < a; < 0.82) Q indicate that 
charge ordering is accompanied by a structural transi- 
tion and antiferromagnetic long range order indeed de- 
velops at lower temperature. The nature of spin fluctu- 
ations changes from ferromagnetic to antiferromagnetic 
at the charge ordering transition B]. Optical reflectivity 
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studies of the same compound for Tco > T > T/v have re- 
vealed the coexistence of a polaron response and a charge- 
gap-Hke structure in the optical response This two 
phases behavior is characterized by domains of ferromag- 
netic and antiferromagnetic spin correlations. Magneti- 
zation and resistivity measurements in Cai_2;Sm2;Mn03 
for < a:: < 0.12 are interpreted in terms of the exis- 
tence of a cluster glass metallic state below some critical 
temperature and demonstrate the lack of true ferromag- 
netism in these electron doped manganites 

From the theoretical point of view, the pioneering work 
of de Gennes proposed a canted phase to resolve the 
competition between the ferromagnetic double-exchange 
(DE) interaction introduced by the presence of itiner- 
ant holes and the super-exchange (SE) interaction. Re- 
cently, several contributions to this problem have been 
reported. Arovas and Guinea (2l], studied this problem 
using a Schwinger boson formalism to obtain a phase dia- 
gram showing several homogeneous phases and pointing 
out that phase separation replaces the canted phase in 
a large region. Indeed phase separation appears in sev- 
eral numerical treatments of the problem p^ . In other 
analytical treatments, more adequate to treat local in- 
stabilities, non-saturated local magnetization states have 
appeared at zero temperature M. Yu. Kagan et. 

al. have studied the stability of the canted phases 
against the formation of large ferromagnetic 'droplets' 
containing several particles and they conclude that the 
formation of droplets is favored in the ground state. The 
variety of results obtained from the different approaches 
points to the need of clarifying the picture and testing 
the results. 

The diversity of results in the phase diagram is con- 
nected to the two dominant magnetic interactions acting 
in these systems: double exchange arising from the Cg or- 
bitals and antiferromagnetic superexchange due mainly 
to half-filled tg orbitals. These competing interactions 
could give rise to textures of different scales according to 
the relation between them. 

Here we study a model for diluted compounds that char- 
acterizes them in terms of a single parameter: the rela- 
tion between superexchange and hopping energies. In 
fact the ionic radii of the intervening dopants varies sub- 
stantially the hopping matrix element through the vari- 
ation of the Mn-O-Mn angle as discussed in Ref. ||3^ . 

In this work, we find the low energy quasiparticles and 
characterize their structure and dispersion relation in the 
low concentration limit. These quasiparticles correspond 
to the electron followed by a ferromagnetic local dis- 
tortion (ferromagnetic polaron) in the antiferromagnetic 
(AF) background. The dispersion relation is dominated 
hy k ^ k + TT scattering due to the presence of AF order. 
In order to make a connection with transport properties, 
we study the tendency to localization of these polarons 
in the presence of impurities and magnetic field. 

We also study the interaction between the quasipar- 



ticles. For two particles, the profile of each quasiparti- 
cle is practically the same as the one obtained for only 
one particle. The effective magnetic interaction between 
these quasiparticles is antiferromagnetic for large values 
of K and becomes ferromagnetic when the size of the po- 
larons increase over some critical value. This change in 
the spin-spin correlation is followed by a change in the 
charge-charge interaction. If the size of the ferromagnetic 
distortion induced by one particle is larger than some 
critical value the other particle could share the same dis- 
tortion giving rise to a bipolaronic bounded state. 

We include diagonal disorder in the two particles prob- 
lem in order to determine how this affects the effective 
interactions between the polarons. We find that ran- 
domness in diagonal energies induce a distribution of ef- 
fective interactions between localized polarons which can 
be ferro or antiferromagnetic. 



II. THE MODEL 

To render evident the nature of the ground state, we 
resort to the Lanczos method, which is free from approxi- 
mations. The Hamiltonian is simplified to a single orbital 
per site, no lattice effects are considered, and we have to 
reduce to one dimensional chains. However our results 
provide a simple picture that, we presume, can put to 
test the dilute limit of electron doped systems. In these 
systems, the limitations of the model Hamiltonian may 
not be as stringent as in the hole doped systems for the 
following reasons: a) the lattice structure is more sym- 
metric so Jahn Teller distortions should play a less im- 
portant role (this observation is consistent with the fact 
that hole doped manganites Lai_a;Caa;Mn03 are insulat- 
ing for low concentrations of dopant while electron doped 
ones Cai_j^Smj;Mn03 are semi-metallic for similar values 



of X 19|); b) the large in-site Coulomb repulsion inhibits 
double occupation so that it may be possible to describe 
the physics by the use of a single effective orbital; c) the 
antiferromagnetic structure of two inter-penetrating lat- 
tices can be properly described in one dimension. 

In order to describe the manganites we consider two 
degrees of freedom: localized spins that represent the t2g 
electrons at the Mn sites, and itinerant electrons that 
hop from Mn orbitals to nearest neighbor e2g or- 
bitals. The model Hamiltonian includes Hund coupling 
( J/i) between localized and itinerant electrons, an antifer- 
romagnetic superexchange interaction between localized 
spins (K), a hopping term of strength t which we will use 
as energy unit hereafter and an on site (U) and nearest 
neighbor (V) Coulomb repulsions: 
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where n^^o- = c+ da {rn = ruT^+nii), and c+ , Cia creates 
and destroys an itinerant electron with spin a at site i, 
respectively. Si and ai are the localized and itinerant 
spin 1/2 operators at site i, respectively. In order to re- 
duce the Hilbert space we take S =1/2 for the localized 
spins instead of S =3/2. It has been shown numerically 
that, in the absence of antiferromagnetic coupling, the 
results for S =3/2 and 1/2 are qualitatively similar [ p2[ . 
In the absence of AF coupling the Hamiltonian reduces 
to the model called Ferromagnetic Kondo Lattice (FKL) 
studied with different numerical methods . The reported 
phase diagram is very similar for dimensions D=l, 2, 3 
and even in infinite dimension showing a ferromagnetic 
phase for large Jh/t |2^i|2^- In the low electron concen- 
tration limit (a; ~ 1) this ferromagnetic phase is obtained 
for any value of Jh- In this limit the AF super-exchange 
coupling between localized spins [K) is needed to get 
the observed AF phase. This model has been studied 
with Monte Carlo method for classical spins and finite 
concentration |^ . Recently some results have also been 
obtained for a very small value of K/t (Kjt ~ 0.05) [p6| . 
In this paper we focus on the dilute limit considering the 
cases of one and two added electrons. Preliminary results 
have been reported in reference |^ . 

In the following we will concentrate on the case U ~ 
y = 0, unless stated otherwise. For the manganites, it is 
believed that the coupling is large. However there are 
some observations in LaMnOs and CaMnOa which indi- 
cate that Jh could be not so large. Quoted values refer to 
Jh ~ ley from optical conductivity data and band- 
width W — leV 1^. In our one dimensional counterpart 
model this would mean Jh = 4t. Band structure calcula- 
tions 1^ also indicate Jh — leV. A larger value of Jh is 
inferred from scanning tunneling spectroscopy pO{ . Re- 
cent optical conductivity for CaMn03 gives 1.7 eV 
and 2eV [Q. However W is also larger in the ferromag- 
netic. Sarma et. al. (3|] indicate W = AeV {t = O.SeV), 
giving Jh/W — 1/2. In one dimension this would corre- 
spond to Jh — 2t. In all cases U = 8 — lOeV » Jh- For 
that reasons we will keep the model general considering 
both the strong coupling {Jh ~ lOeF) and weak coupling 
{Jh ^ leV) regimes. 



the particle can be determined from correlation functions 
which mix charge and spin variables. In this way it is pos- 
sible to distinguish between polaronic and non-polaronic 
regimes. 

The quasiparticle character of this polaronic distortion 
becomes evident from the dispersion relation. The effec- 
tive mass of the quasiparticle is inversely proportional 
to the bandwidth. In this way it is possible to study 
the variation of the effective mass as a function of the 
parameters. 

By lowering the diagonal energy in one site it is possible 
to localize the polaron. We study the dependence of the 
localization length and the spin distortion on the effective 
mass and on an external magnetic field. 

The correlations functions of this section have been 
calculated using periodic boundary conditions (PBC). 

A. Polaron profile 

In order to determine the polaron profile we calcu- 
late < UiSjSj+i > for the ground state. Because of 
translation symmetry this correlation function depends 
only on |i — j|. The results for different values of K/t 
and Jh — lot are shown in Fig.|l|(a) where we plot 
N < TiQSjSj+i > vs j, where N is the number of sites. 
As it can be seen in Fig.|l|, for large j this correlation 
function takes a value very close to the one obtained 
from the Bethe ansatz solution of the Heisenberg chain, 
< SjSj+i >= —0.443. The extension and the magnitude 
of the spin distortion around the particle increases as K/t 
decreases. For large values of K/t, the background of lo- 
calized spins remains in a 5 = total spin state although 
it presents distortion in the spin-spin correlation function 
due to the itinerant particle. By decreasing K/t the total 
spin increases to 5 = 1 (this occurs for K/t w 0.3 in the 
case of |l|(a) for Jh/t = 10) and higher values indicating 
a net polarization of the localized spins. In this case, we 
can refer to a polaronic regime. 



III. ONE PARTICLE 



In this section we investigate the homogeneity of the 
solutions for different values of K/t. To this end we cal- 
culate the ground state with one itinerant e2g electron 
added, for chains of different sizes up to TV = 20. For 
large values of Jh , the particle modifies substantially the 
spin structure in its surroundings forming magnetic po- 
larons. The distortions of the magnetic structure around 
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FIG. 1. a) We show the correlation function < noSjSj+i > 
for a 16 sites chain, with = 10 and different values of K 
{K = 0.1 triangles, K — 0.2 open squares and K = 0.3 full 
squares) The maximum value of 5 for each value of K 

is indicated in the figure. One can observe that the total 
S and the extension of the magnetic distortion increase as 
K decreases. The oscillatory behavior found at if = 1 is a 
consequence of the weakening of the antiferromagnetic links 
around the charge. The inset shows < SjSj+i > for a 16 
sites Heisenberg chain where the link at site zero is a factor 2 
smaller than the rest. 6) The same as a) for the weak coupling 
regime (A = 1). 

In Fig j|(b) , we show the same correlation function for 
Jh ~ t. For this value of Jh we are already in the weak 
coupling hniit where the spin distortion is almost null for 
any value of K/t (notice the change in the vertical axis 
scale). The profile of the distortion in this case is, in fact, 
very similar to the one obtained in the strong coupling 
limit {Jh ~ lot) for large values of K/t. The polaronic 
regime appears, in this case, for much lower values of 
K/t. 

The oscillations observed in the curve corresponding to 
K = 0.3 are also observed for larger values of K. They 
are a consequence of the weakening of the antiferromag- 
netic links around the charge position which produce a 
sort of local spin dimerization. Because of the compe- 
tence between double-exchange and super-exchange, the 
effective antiferromagnetic interaction between the site 
where the particle is and the nearest neighbor one is 
weakened {{noSoSi) > —0.443). The next spin corre- 
lation {{nQSiS2)) may be even stronger than in the un- 
doped case {{noSiS-z) < —0.443) due to this weaken- 



ing. This forces the next link ((no<5'2S'3) > -0.443) to 
be weaker, and the same kind of reasoning can be ap- 
plied to the rest. This effect explains the oscillations of 
(uQSjSj+i) as a function of j for large values of K/t. 
To prove this point we show in the inset the nearest 
neighbors spin-spin correlation functions for a Heisen- 
berg chain of the same size where the link between site 
zero and one is a factor of two smaller than the rest. 
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FIG. 2. a) Evolution of the polaron profile from the weak 
coupling regime {Jh = 1) to the strong coupling limit 
{Jh = 10). b) Magnetic correlation between the site where 
the particle is and the nearest neighbor as function of A. 

In Fig.^(fe) wc show the variation of the polaron profile 
as a function of Jh, for K/t = 0.1. In this way it is 
possible to see the transition from non-polaronic (weak 
coupling) to polaronic (strong coupling ) regime. 



B. Dispersion relation 

As K decreases, it is difficult to find an adequate ap- 
proximation to describe the large polaronic distortion. In 
order to obtain the effective mass of these polarons, we 
investigate the dispersion relation for charge excitations. 
To this end we calculate the lowest energy state for dif- 
ferent values of the momentum k = 2T:n/N within the 
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subspace where the total spin is that of the ground state. 
In Fig. 2 from Ref. ||2^] we show the dispersion relation 
scaled to the thermodynamic limit for K — 3, Jh — 100; 
K = 1, Jh = 10; and K ^ 0.3, A = 10. 

We start analyzing the dynamics in the regime where 
(J/i >> K >> t). In this case, the charge moves as a 
spin one (S). The effective hopping resulting from the 
projection of the hopping term onto the reduced S — 1 
Hilbert space is: tPijCEiSj + 1/2), where Pij is the per- 
mutation operator between sites i and j. We can picture 
the movement of the particle, in this limit, as going from 
a state iTitfititoan intermediate state i T i T 
t i) and finally to i t i T I tT i, where (0) repre- 
sents the Sz = +1(0) components of the spin S = 1. 
Thus, in order to move, the charge has to hop to the 
nearest neighbor, via a spin flip process, through states 
that differ in energy by A « K/2. It can be easily veri- 
fied that the effective hopping of this process is equal to 
tef = t/^/2. The dispersion relation given by this dynam- 
ics is: A/2 ± ^(A/2)2 -|-4i^^cos2(fc) . The expression 

corresponding to the lower band fits well the numerical 
results ioi K = 3 and Jh — 100 as shown in Fig. 2 from 
Ref. . This expression is valid in general for a particle 
moving in an antiferromagnetic background where scat- 
tering between k and k+ir states dominates the dynamics 
of the particle. 



T ' 1 ' r 




1/N 

FIG. 3. Scaling of the difference between fc = and k = n 
states for different values of K and Jh ~ 10. 

In Fig.^ we show the scaling of the energy difference 
between k and k + n states. This difference must tend to 
zero as N tends to infinity due to the reasons explained 
above. In this way we can check the quality of the scaling 
used to calculate the points in Fig. 2 from Ref. The 
linear extrapolation goes to zero with an error lower than 
5 X 10^"^ for the three values of K/t considered. 

In the case K » Jh » t the spin distortion can be 



neglected and the particle propagates in an antiferromag- 
netic lattice. The Hund interaction alternates the site en- 
ergy of the propagating particle so that the difference be- 
tween the two sublattices is given by A = Jh{< (TjSj+i > 
— < ajSj >) = Jh{< SjSj+i > — < OjSj >) where we 
approximate < (JjSj >~ 1/4 its value at the triplet state, 
and < cTjSj+i >« < SjSj+i >= In 2 - 1/4, the Bethe 
ansatz value. Using these values we find A = 0.19 J/,.. In 
this case t^f is equal to t. 

When Jfi >> t > K, the magnetic distortion around 
the charge is large and the effective hopping is domi- 
nated by the overlap between the magnetic distortions 
about the nearest neighbors sites. This last effect dom- 
inates the polaron effective mass. Therefore, the mass 
of polarons increases when K decreases, as in Fig. 2 from 
p7| , where tef decreases from 0.75 for K = 1 to 0.23 
for K = 0.3, showing that the spin distortion around 
the charge increases in magnitude and extension when 
K decreases. 

The bandwidth calculated for several values of K (Fig. 
3 of Ref. H^]) clearly shows two regimes: K « t and 
K » t. The first corresponds to a large magnetic dis- 
tortion and the second corresponds to a smaller one ac- 
cording to Fig.|l|. The bandwidth goes to zero with K in- 
dicating that the effective mass of the polaron increases 
continuously by decreasing K . When K is equal to zero 
the size of the ferromagnetic distortion becomes infinite 
and we can not consider any more this distortion as part 
of the quasiparticle structure. The quasiparticle scenario 
is no longer valid (it is not possible to associate a mag- 
netic distortion to the electron since the magnetic dis- 
tortion is infinite). In this limit, the electron moves in 
a ferromagnetic background freely and thus generates a 
bandwidth equal to the tight-binding At. In our case, 
since we have considered finite systems, this happens for 
a small, but not zero, value of K . This regime was not 
shown in Fig. 3 of Ref. . 

C. Localization 

To test how robust the polaronic description is, we 
pin the polaron by changing in —eg the diagonal energy 
at site zero. This may be relevant to the real materi- 
als since the doping process necessarily introduces some 
disorder. This always localizes the particle in a linear 
chain, but the localization length should be very differ- 
ent for different effective masses. A small eo localizes 
much more the polaron for low values of K than for larger 
ones. This is shown in Fig.^(a) where we plot < rii > 
around site zero for different values of K, while fixing 
Jh = lot and eo = 0.05. We also show in Fig.|(b) the 
correlation function < noSjSj^i > in order to determine 
the magnetic distortion associated with each localized 
state. While charge localization does not change very 
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much for K between 0.02 and 0.1, the magnetic distor- 
tion increases monotonicaUy by decreasing K. To under- 
stand this point is useful to take the hmit K/t — and 
then slowly increment K/t. For K/t = 0, the system is 
a fully polarized ferromagnet and thus the ferromagnetic 
distortion is infinite. In this case, however, the localiza- 
tion length is finite, although the largest possible. As 
K/t increases the magnetic distortion will shrink, with- 
out affecting the localization length, at least until they 
are comparable. In this regime the magnetic distortion is 
decreasing while the localization length is not changing. 
The difference between the localization and the magnetic 
distortion lengths depends on the value of K and eo, and 
determines the limit of validity of the quasiparticle con- 
cept. 
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FIG. 4. Charge localization, a) We show < > as a 
function of position in a 16 sites chain where the energy at 
site zero is changed from the rest by — eo=-0.05 units of t. 
The short localization length of the lower K curves indicate 
effective bandwidths of the order of the energy change. 6) We 
show the extension of the magnetic distortion induced by the 
localized particle. 



localization length starts to decrease as a consequence of 
the increase of the effective mass. This behavior remains 
up to some value of K where the size of the magnetic 
distortion is similar to the localization length. If we fur- 
ther decrease the value of K then the localization length 
starts to increase up to the value corresponding to the 
free case for K << eq, and we loose the quasiparticle 
picture. 
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FIG. 5. Same as Fig.| 
(A = 1). 



for the weak coupling regime 



In Figs|5|(a) and (6) we show the same as in Fig.^ but 
for the weak coupling regime (J/i ^ t). In this case the 
localization length is always larger than in the strong 
coupling regime as expected. The size of the magnetic 
distortion is also smaller. 



For large values of K the polaronic distortion is very 
small and the localization length does not change sub- 
stantially because the effective mass of the quasiparticle 
is similar to the free mass. When K becomes smaller the 
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FIG. 6. Effect of magnetic field on diarge localization in 
the strong coupling limit {Jh = 10). We show a) < rii > 
for different values of an external magnetic field for the same 
chain as in Figj^; b) the evolution of the magnetic distortion 
as function of the magnetic field. 

It is interesting to note the differences in the profile 
of the magnetic distortion between the localized case 
(Fig.^(a)) and the homogeneous case (Fig.|l|(a)). In the 
localized case the magnetic distortion is abrupt defining 
the polaron limits quite sharply (one could think of a 
step shape). In the non- localized case the polaron profile 
decays more continuously defocusing the polaron limits 
(one could think of an exponential- like shape). 

In Fig 8(a), we show the change in the values of < > 
for different magnetic fields in the strong coupling regime 
{Jh = lOi and K = O.lt). Sz increases by one between 
two successive values of H starting from Sz — 2.5 for 
H = 0. It can be seen that the localization of the po- 
laron decreases with magnetic field as a consequence of 
increased effective hopping between nearest neighbors. A 
fact that may be important for the transport properties of 
these systems since it implies a negative magnetoresistive 
behavior for conductivity due to hopping between local- 
ized states |3^Jl^,|l^ . The size of the magnetic distortion 
relative to the background increases with the magnetic 
field (see Fig|6|(b)). This could be due to the stronger 
effect of the magnetic field on the spins which are close 



to the boundary of the polaron. Those spins are form- 
ing weak links due to the competition between super and 
double exchange so they must have a larger susceptibility. 
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FIG. 7. The same as Fig.^ in the weak coupling limit 
(A = 1). 

In Figj^(a) we show the effect of magnetic field on the 
localization for the weak coupling limit (Jh — 1, K = 
0.1). It is clear from this figure that the magnetoresistive 
effect is very weak compared with the strong coupling 
limit. In addition, the size of the magnetic distortion is 
very small (Fig.0(6)) as expected for the weak coupling 
regime. 



IV. TWO PARTICLES 

In this section we study the interaction between the 
quasiparticles. To this end we calculate the magnetic dis- 
tortion induced by two itinerant electrons, as well as the 
binding energy and the charge-charge correlation func- 
tion. In this way it is possible to discriminate if the 
effective interaction is attractive or repulsive. 



7 



LA-UR-99-3625 



August 24, 2000 



accepted in Physical Review B 




J 

FIG. 8. Charge-charge correlation as a function of distance 
in the strong coupling limit (Jh = 10). 

To obtain closed shell conditions in the spinless limit 
{K/t = and Jh — oo) we have taken antiperiodic 
boundary conditions (APB). We have only used PBC for 
the calculation of binding energies because this quantity 
also involves the energy of one particle added, which has 
been calculated with PBC. 

The charge-charge correlation function, (norij) gives 
information about the character of the interaction be- 
tween the polarons. This correlation function is shown 
in Fig.^ for a chain of 14 sites. The maximum always 
occurs for the largest possible distance between the par- 
ticles, what could indicate a tendency towards repulsion. 
However in some cases this could just indicate that the 
size of the chain is similar or smaller than the mean sep- 
aration between the two particles in a bound state. To 
discriminate between these two possibilities it is neces- 
sary to calculate the binding energy for different sizes of 
chains. 
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FIG. 9. Binding energy, (£(2) + £^(0) - 2£(1)), as function 
of the inverse of the system size in the strong coupling limit 
{Jh = 10). 

Fig.^(a) shows the binding energy {E{2) + E{0) — 
2E{1)) as a function of the inverse of the chain length, 
in the strong coupling limit. We find that, although up 
to the sizes we can compute the interaction is always 
repulsive for K > 0.1, the extrapolation to the thermo- 
dynamic limit seems to give a negative value, suggesting 
the possible existence of a bound state. When K is low 
enough (K = 0.05t) the binding energy is already nega- 
tive for > 12. The fact that the interaction is repulsive 
for chains smaller than 12 sites indicates that the mean 
separation of the bound pair is around six lattice param- 
eters. In Fig.^(b) we show the weak coupling limit. Now 
for every value of K we get a vanishing binding energy 
extrapolated to the thermodynamic limit. This reflects 
the importance of the double exchange in the dynamics 
of the electrons. By increasing the coupling the system 
seems to cross between a non-bounding regime into an 
attractive one. 



8 



LA-UR-99-3625 



August 24, 2000 



accepted in Physical Review B 





□ — □ 



o— o 



4 



—r- 
6 




J 

FIG. 10. Profile of the the magnetic distortion with two 
particles added for different values of K. 

In Fig.^we show the profile of the magnetic distortion 
for different values of K and for the maximum separation 
between the particles. As for the one particle case, the 
size and amplitude of the magnetic distortion increases 
as K decrease giving rise to a strong polaronic regime 
for K ^ O.lt — 0.2t. The two polarons are practically in- 
dependent of each other, since their shape is essentially 
the same as in the single polaron regime described in 
the previous section (see Fig.Q). If the size of each po- 
laron becomes larger than some critical length, one would 
expect that the second particle added could gain some 
magnetic energy by sharing the distortion created by the 
first one. This mechanism is similar to the spin-bag idea 
proposed by Schrieffer to explain the attractive in- 
teraction giving rise to pairing in the high temperature 
superconductors. From this consideration, the existence 
of a bound state under some critical value of K can be 
understood. 




FIG. 11. Spin-spin correlation as function of distance for 
the largest separation between the charges. 

If the spin-bag picture is correct, the polarons must in- 
teract ferromagnetically when they form a bound state. 
This can be checked by calculating the total spin of the 
ground state and the spin-spin correlation between the 
charges. In Fig, 11 we show spin-spin correlation as a 
function of distance when one particle is at site zero and 
the other one is at the longest distance in the chain (site 
7 in a 14 sites chain). It is clear that, as K decreases, 
the interaction between the two charges becomes more 
ferromagnetic. To confirm this result, in Fig.|l2|, we show 
the spin-spin correlation function between the two par- 
ticles from where the effective magnetic interaction be- 
tween the two quasiparticles can be extracted. It is clear 
that this interaction evolves from a very small value for 
K — 0.4 to the largest ferromagnetic value for K — 0.1. 
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FIG. 12. Spin-spin correlation between the two added elec- 
trons for different values of K. 

This change in the behavior of the interaction is also 
obtained in the charge-charge correlation function (Fig.|^) 
where, for some critical value oi K between 0.1 and 0.05t, 
we observe a change in the variation of the intensity of 
the maximum with K. It increases with decreasing K 
from 0.3 to 0.1 and then it decreases practically the same 
magnitude between K = 0.1 and K = 0.05. 




J 

FIG. 13. Spin-spin correlation function between the two 
particles as function of the relative distance. The diagonal 
energies of sites and 7 have been lowered by eo in order 
to pin the two quasiparticles for K = 0.2. eo=0.01 (solid 
squares), eo = 0.1 (open circles), eo = 1 (solid circles). 

We also studied the interaction between polarons in 
presence of disorder. As it is seen in Fig.^, lowering of 
the diagonal energy of one site not only pins the quasi- 
particle, but also changes its size. This change should in- 
duce a modification of the effective magnetic interaction 
between two pinned quasiparticles. This can be seen in 
Fig.|l^ where we show the spin-spin correlation between 
the two particles a function of distance and for different 
values of cq. The diagonal energy has been changed by 
— eo in both sites and 7. The effective magnetic inter- 
action evolves from ferromagnetic values for eg < 0.1 to 
antiferromagnetic for eo > 0.1. This can be easily under- 
stood from the reduction in the size of the polaron dis- 
tortion which is induced by the increment of the pinning 
energy eo. As each particle is more localized, the over- 
lap between polaronic distortions decrease. Each particle 
cannot see the magnetic distortion generated by the other 
and the effective magnetic interaction is dominated by 
the superexchange through the localized spins between 
the two polarons. 




J 

FIG. 14. Same as Fig.|l^, but for different values of K and 
eo — 0.05t. K — 0.05 (solid circles), K = 0.1 (open squares), 
K = 0.2 (open circles), K — 0.3 (solid triangles), and K = 0.5 
(solid squares). 

In Fig.0 we show the variation of the effective mag- 
netic interaction between two pinned quasiparticles as a 
function of K. In this case we get a non-monotonic be- 
havior. By decreasing K from 0.5 to 0.2 the effective in- 
teraction becomes more ferromagnetic. But for K = 0.1 
this behavior changes drastically becoming antiferromag- 
netic when the particles are close to each other and fer- 
romagnetic when each one is close to its pinning center. 
This indicates the formation of two well defined ferromag- 
netic distortions surrounding each center of pining with 
an antiferromagnetic effective interaction. For K = 0.05 
this effective interaction becomes ferromagnetic and we 
get a uniform ferromagnetic state (the size of the induced 
distortion is equal to the size of the system). 

Localizing the particle reduces the overlap between po- 
larons, the system becomes more inhomogeneous and the 
effective interaction between the polarons evolves from 
ferromagnetic to antiferromagnetic. For example, for 
K = 0.1, (SoSjUorij) is positive with eo = (FigJTl|) 
and becomes negative for eo = 0.05 (Figjl^). 

Finally we have studied how sensitive are the results to 
the inclusion of U and V in the strong coupling regime. 
We do not find any significative change in the physics due 
to the presence of these interactions. As we are study- 
ing the dilute limit, this interactions should not be so 
important as for higher concentrations. Even though, if 
both particles are at the same site, they must be in a 
singlet state due to the Pauli principle. This state has 
a very large energy (~ Jh) in the strong coupling limit 
{Jh = 10) giving rise to a very small double occupation 
in the ground state. This explains the poor role of U. To 
explain the weak effect of V it is useful to take a look at 
Fig.^ where it can be seen that the occupation of adjacent 
sites for the two particles is very small. 
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V. CONCLUSIONS 

In summary, we have investigated the possibihty of 
strong magnetic distortions induced by carriers in the 
ground state of the DE-SE model Hamiltonian. The 
model describes chains of localized spins coupled anti- 
ferromagnetically and interacting ferromagnetically with 
itinerant electrons. We have studied in detail the case 
with one and two electrons. Assuming the model ade- 
quately describes the physics of electron doped mangan- 
ites, the results presented here point to a picture of these 
systems where heavy polarons dominate the magnetic 
and transport properties. Their masses depend strongly 
on the relation between the hopping energy and the an- 
tiferromagnetic superexchange interaction. Clearly, the 
doping itself will localize the polarons so that transport 
will result from hopping between pinned sites. Negative 
magnetoresistance should appear as a consequence of the 
decrease of the pinning energy with magnetic field |3^ ] . 

We have done scaling on the binding energy of two 
polarons. Two regimes can be identified depending on 
the coupling (A). In the strong couphng case, although 
only for small values of K we find a negative binding 
energy {K < 0.05 and N > 12), the scaling to the ther- 
modynamic limit is negative for every K we have studied, 
indicating the possible existence of a bound state. In the 
weak coupling case, for most values of K, the scaling to 
the thermodynamic limit is zero, and all values are posi- 
tive. However, as the charge-charge correlation function 
(see Fig.||) is always peaked at the largest possible dis- 
tance between the particles, it is necessary to confirm our 
scaling results in larger systems. 

The existence of a bound state can be understood from 
the ferromagnetic distortion of the polarons. Below some 
critical value of K this ferromagnetic distortions of two 
neighbor polarons overlap and the charge-charge inter- 
action becomes attractive. This scheme resembles the 
spin-bag picture proposed by Schrieffer to explain the 
pairing mechanism in the High Tc compounds. The effec- 
tive magnetic interaction between polarons also changes 
from antiferromagnetic to ferromagnetic by decreasing 
the value of K. 

We have also studied the effect of localization on the 
effective interaction between the quasiparticles. The ef- 
fect of localization is not only to pin the quasiparticles 
but also to change the shape of the ferromagnetic distor- 
tion. As the shape and size of this distortion is clearly 
connected with the spin-spin and charge-charge effective 
interaction we get that the amount of disorder play a very 
relevant role in the low energy physics of these quasipar- 
ticles. Depending on the pining energies, ferromagnetic 
and antiferromagnetic interactions could be generated. 
In real systems frustration induced by such competing 
interactions could give rise to spin glass behavior. In 
this way, the interplay between diagonal disorder and 



effective magnetic interactions between polarons could 
explain the cluster glass character of the metallic phase 
found by Maignan et. al. [|l^ in Cai_2;Sm2;Mn03 for 
< a; < 0.12 . 

Finally, we would like to point out that the order of 
oxygen vacancies in CaMnOa-^ makes real the possibil- 
ity of one dimensional electron paths in these materials 
p6[ . We hope that our results will stimulate more ex- 
perimental and theoretical investigations on the electron 
doped manganites. 
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